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There are two distinct techniques of proven effectiveness for extending the coherence lifetime of 
spin qubits in environments of other spins. One is dynamical decoupling, whereby the qubit is 
subjected to a carefully timed sequence of control pulses; the other is tuning the qubit towards 
‘optimal working points’ (OWPs), which are sweet-spots for reduced decoherence in magnetic fields. 
By means of quantum many-body calculations, we investigate the effects of dynamical decoupling 
pulse sequences far from and near OWPs for a central donor qubit subject to decoherence from a 
nuclear spin bath. Key to understanding the behavior is to analyse the degree of suppression of the 
usually dominant contribution from independent pairs of flip-flopping spins within the many-body 
quantum bath. We find that to simulate recently measured Hahn echo decays at OWPs (lowest- 
order dynamical decoupling), one must consider clusters of three interacting spins, since independent 
pairs do not even give finite T 2 decay times. We show that while operating near OWPs, dynamical 
decoupling sequences require hundreds of pulses for a single order of magnitude enhancement of T 2 , 
in contrast to regimes far from OWPs, where only about ten pulses are required. 

PACS numbers: 03.67.Lx,03.65.Yz,76.60.Lz,76.30.—v 


I. INTRODUCTION 

Individual electronic and nuclear spins are promising 
candidates for realizing scalable quantum computing in 
solid state systems such as silicon.® In addition, they of¬ 
fer a valuable test-bed for the experimental investiga¬ 
tion of decoherence driven by quantum baths which typ¬ 
ically comprise large numbers of spins surrounding a cen¬ 
tral spin qubit. Despite the large number of bath spins 
involved, accurate simulation of experimental coherence 
decays is computationally tractable since the bath can of¬ 
ten be decomposed into independent contributions from 
many small sets or clusters of spins.®® In many cases, 
the dominant contribution to the decoherence dynam¬ 
ics arises from pairs of bath spins; in effect, from the 
magnetic noise due to the independent ‘flip-flopping’ of 
spin pairs. Contributions from larger clusters are usu¬ 
ally only needed for high accuracyp*® however, there is 
also interest in cases where experimental evidence for the 
many-body or large-cluster nature of the bath might be 
most evident. For example, it was recently found that 
applying dynamical decoupling sequences of microwave 
pulses may enhance many-body correlations, relative to 
the independent spin-pair contribution.® 

A defining characteristic of a quantum bath is the 
back-action between the central system and environ¬ 
ment, making the bath dynamics sensitive to the state 
of the central spin. A particularly striking example 
has been identified recently in spin systems with special 
coherence sweet-spots termed ‘optimal working points’ 
(OWPs),fi2HIll where theory and experiment found that 
coherence times can change by orders of magnitude with 
even small variations in applied magnetic field 200 
G). A drastic change in back-action occurs via changes to 


the states of the central spin only: the change in external 
field has little direct effect on the bath spin-pair dynam¬ 
ics. However, it was found that usual cluster simulations 
using only pairs predicted infinite coherence times for the 
lowest order dynamical decoupling sequence (Hahn spin 
echo), at odds with experiment.^ 

These findings motivate us to investigate here the in¬ 
terplay between dynamical decoupling and OWPs for a 
quantum bath system and in particular, to clarify where 
and to what extent, the independent spin-pair contri¬ 
bution dominates. Dynamical decoupling is one of the 
most established methods for extending coherence.^i^^f^ 
It involves subjecting the qubit spin to a sequence of mi¬ 
crowave or radio pulses. A wide variety of solid state spin 
qubits have been studied under dynamical dec oupling 
control; these include Group V don ors i n silicon,®2SH2ll 
nitrogen vacancy centres in diamond,®^^!^ GaAs quan¬ 
tum dote^ rare-earth dopants in silicates ,l^malonic acid 
crystal^^ and adamantane.®^ 

Long-lived coherence is a key requirement for imple¬ 
menting fault-tolerant quantum computation,^^ as well 
as quantum memory.l^^ One way of extending coherence 
times T 2 in silicon is to use isotopically enriched samples 
in which the abu ndance of nuclear spin isotopes is sig¬ 
nificantly reduced.However, it is advantageous to 
retain nuclear spins for their potential use as long-lived 
quantum registers.^ It is thus also of practical impor¬ 
tance to understand whether dynamical decoupling and 
OWP techniques may be advantageously combined for 
a quantum bath of nuclear spins. For donor electronic 
qubits in silicon, it is known that due to inhomogeneous 
broadening from naturally-occurring ^®Si spin isotopes, 
there is a significant gap bet ween the T 2 ^ 100 ms in 
natural silicon near an OWP,li®i^ and the T 2 ~ 2 s in 
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isotopically enriched ^®Si with a low donor concentration 
at the same OWP.^ Also, dynamical decoupling may be 
useful when it is convenient to operate with the magnetic 
field close to but not exactly at the OWP. 

It is well established that for dynamical decoupling to 
be effective, the pulse spacing r = t/2N for a sequence 
of N control pulses (where t is the total evolution time) 
cannot exceed the correlation time of the bath noise. But 
the relevant correlation time, in turn, is an emergent 
property of the underlying microscopic quantum bath, 
comprised of typically ^ 10"^ — 10^ significant clusters of 
spins of different coupling strengths, different sizes and 
subject to varying degrees of back-action from the cen¬ 
tral qubit. Therefore, to quantitatively simulate the re¬ 
sponse to dynamical decoupling, a realistic simulation of 
the combined system-bath dynamics at the microscopic 
level is important. 

In this paper, we present quantum many-body sim¬ 
ulations of the system-bath dyn amics using the cluster 
correlation expansion (CCE),SE5I including contributions 
from clusters of up to 5 spins (CCE5) as illustrated in 
Fig. [3b). We compare coherence decays at an OWP 
with regimes far from an OWP (denoted by ‘y^OWP’) 
with a view to identify the optimal strategy for enhancing 
electronic spin coherence times of donor spins in natural 
silicon. Simple analytical expressions for the behavior of 
independent bath pairs coupled to the qubit aid under¬ 
standing in all the regimes we consider. 

In addition, our work fills an outstanding gap in un¬ 
derstanding decoherence near OWPs. It was recently re¬ 
ported that near OWPs, numerical calculations of the 
Hahn spin echo {N = I) involving only independent pair 
dynamics (CCE2) yield results in conflict with experi¬ 
ments: coherences decay initially, then after a short time, 
the decays stop.ff^In the present work, we undertake the 
computationally more challenging CCE3 — 5 many-body 
calculations in order to clarify the origin of the measured 
coherence decays. We find that three-spin clusters suf¬ 
fice to give decays in good agreement with previously 
reported experimental results.^ 

It is worth clarifying the physical meaning of the 
above-mentioned three-cluster result. It is not a matter 
of enlarging the quantum bath with additional nuclear 
spin clusters of the same size. As illustrated in Fig. 0b). 
a three-spin cluster (blue) can be decomposed into three 
distinct flip-flopping pairs (each nuclear spin can con¬ 
tribute to more than one flip-flopping pair). Put simply, 
if all such three-clusters in a given, randomly generated 
set of impurities in a crystal are decomposed into the 
constituent flip-flopping pairs, an infinite decay time is 
obtained. If, however, the exact same configuration of 
spin impurities are aggregated into the ‘triangle’ struc¬ 
tures illustrated in Fig. int)- the correct experimental 
behavior emerges. To our knowledge, there is no other 
example of a central spin system which so fully eliminates 
the pair-driven dynamics. 

In contrast, for t^OWP regimes, Hahn decays for 
dipole-allowed transitions are well described by CCE2.I^ 


However, for modest A^ < 10 pulse numbers there can be 
a large correction to CCE2 from clusters of 3 — 4 spins 
and for even N, but CCE2 still gives T 2 on the correct 
experimental timescale.l^For larger N, we find that CCE 
numerics including only independent pairs (CCE2) once 
again gives converged decays in all regimes whether in 
OWP or T^OWP regimes, so many-body calculations be¬ 
come progressively less important as N ^ 00 . 

The paper is organized as follows. In Section H we re¬ 
view the theory of decoherence and mixed donor qubits 
possessing OWPs. Our main results are presented in Sec¬ 
tion HI, first for the Hahn spin echo, then for dynamical 
decoupling with moderate N. In Section IV, we provide 
an analysis on the suppression of the independent pairs 
contribution to decoherence as V —>■ 00 , as well as the 
magnetic field approaching the OWP: B —)• Howp- Fi¬ 
nally, we present large N results in Section V and discuss 
important aspects of our results in Section VI. 

For our dynamical decoupling calculations, we have 
chosen the Carr-Purcell-Meiboom-Gill (CPMG) se¬ 
quence which applies a set of N periodically spaced 
near-instantaneous pulses (GPMGV) as illustrated in 
Fig.[3d).ESlIZl2ll 

II. DECOHERENCE AS QUBIT-BATH 
ENTANGLEMENT 

A. Overlap of conditional bath evolutions 

For quantum baths, decoherence of a central spin sys¬ 
tem is understood in terms of entanglement between the 
system and the bath, as the combined system evolves 
under a total Hamiltonian given by: 

Htot = HqS “b Hint + .ffbath, (1) 

where Hqs is the central system Hamiltonian including 
all internal nuclear and electronic degrees of freedom, 
while i?bath is the bath Hamiltonian and Hint describes 
the system-bath interaction. 

We start by placing the initial system state in a co¬ 
herent superposition of an upper state \i = u) and a 
lower state \i = 1) by applying a 7r/2 pulse: |4'(t = 0)) = 
(|u) -I- |/)) 0 \B{t = 0)), where \B{t — 0)) is a prod¬ 
uct state of the eigenstates of the non-interacting bath. 
Under the joint system-bath dynamics, the initial state 
evolves into an entangled state: 

V2 |vl/(t)) = \u) 0 \Buit)) + \l) 0 \Biit))) , 

( 2 ) 

where \Bu,i{t)) = Tu,i \B{t = 0)) with unitaries Tu,i for 
each of the two system levels, and we have assumed no 
depolarization of the central states during the evolution. 
The measured temporal coherence decays |/l(t)| can be 
simulated if one can accurately calculate the resultant 
overlap between the bath states correlated with the upper 
and lower qubit states: 

\Cit )\cx \{Buimm = |(S(0)|ftf,|e(0))|. (3) 
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FIG. 1: (color online) (a) The spectra of donor spin systems such as arsenic, antimony or bismuth (pictured) are affected 
by strong mixing between the electron and host nuclear spin, at magnetic fields B smaller or comparable to the hyperfine 
coupling A, allowing a richer behavior than unmixed electron spins. At particular field values termed optimal working points 
(OWPs), decoherence can be strongly suppressed; the arrows indicate the transitions with four of the most significant OWPs. 
(b) Coherences of the central electronic spins are dephased primarily by a surrounding quantum bath of clusters of 2, 3, 4 
or more nuclear spin impurities (for natural silicon, pictured) or other donors (for isotopically enriched silicon), (c) Effect of 
dynamical decoupling (CPMG with even pulse numbers N) as N ^ oo. Plots showing enhancement of the electron 

spin coherence time T 2 as a function of pulse number N, relative to the N = 1 Hahn echo value. We find that while dynamical 
decoupling far from the OWP enhances T 2 by an order of magnitude with about 10 pulses, in contrast, behavior close to an 
OWP is insensitive to dynamical decoupling for low N. For high N, the behaviors near and far from OWPs become comparable, 
(d) Illustrates coherence enhancement as B —> Bqwp (the Hahn spin echo time is plotted). The OWP is for a bismuth 
donor in natural silicon, investigated experimentally in Refs. [T?I13I Inset of (d): The CPMG dynamical decoupling sequence 
consists of the initial 7r/2 pulse, followed by the —t — tv — r—echo sequence repeated N times. The coherence times in (c) are 
when the CPMG decays in Fig. and the fits to the decays in Fig. have fallen to 1/e. OWP results are for the |14) —> |7) 
transition for which Bowp = 799 G. For the field value near the OWP {B = 795 G) in (c), ~ 96 ms while ~ 0.79 ms 

in the t^OWP regime {B = 3200 G). The OWP curve in (d) was calculated using the analytical formula Eq. Q. 


Since the initial bath states are usual trivial thermal spin 
states, the challenge is to obtain the corresponding uni- 
taries Tu^i for extremely large baths (« 10"^ spins for con¬ 
vergent CCE). 

OWPs are sweet-spots in B-field values where the two 
unitaries involving the upper and lower levels equalise: 
T) ~ Tu, occurring when ~ This means that the 
state given by Eq. ([^ can be written as 

|vl/(t)) = |u) + 10) ® T„(t) 1^(0)), (4) 

with the product form preserved (i.e. the state is 
no longer entangled) and therefore decoherence is 
suppressed: |/:(t)| cx \{B{0)\fl{t)fu{t)\B{0))\ = 

|(S(0)|i|S(0))| = l. 


B. Interaction and bath Hamiltonians 

We consider the situation where the central spin inter¬ 
acts with a spin-1/2 bath (e.g. ^®Si impurities) primarily 
through the contact hyperfine interactionP 

ffint JaS-ia, (5) 

a 

where S represents the central electron spin, Ja is the 
strength of the contact hyperfine interaction and a la¬ 
bels the bath spins la- The bath Hamiltonian consists of 


nuclear Zeeman terms and dipolar coupling among bath 
spins: 

^^bath = ^^D+^^NZ) 

Hnz = 

a 

i?D = ;^iaD(ra6)ib, (6) 

a<6 

where ■jn is the gyromagnetic ratio of bath spins and 
D(ra&) is the dipolar tensor coupling bath spins a 
and b separated by Tab, with components Dijir) = 
- SviTj/r^), where i,j = {x,y,z}, 6ij is 
the Kronecker delta and yo = 10“^ NA“^ is the mag¬ 
netic constant divided by AttW^ 


C. Example central system with OWPs 

The dephasing properties of the central spin qubit and 
bath ar e ext remely well studied for the case of a spin-1/2 
qubit ptHHHISI Good agreement with experiment has been 
achieved by cluster-based methods such as the CCE,I^ 
which decompose the bath dynamics into products of 
contributions from clusters of 2, 3, or more interacting 
bath spins as illustrated in Fig. [^b). But simple spin- 
1/2 systems do not have OWPs. 

For the donor systems, however, the central spin 
Hamiltonian Hcs — ujqSz + Alh ■ S, contains the usual 
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Zeeman term arising from the external magnetic field 
(wo = jeB, where 7 e is the electronic gyromagnetic ratio) 
and also a significant hyperfine coupling (A) of the host 
spin to the electron. For example, for the bismuth 
donor, Ih = 9/2 and A = 1.475 GHz, thus the mixing 
between host nuclear and electronic spins becomes sub¬ 
stantial for H ~ 0 — 0.3 T as seen in Fig. [^a). 

Details of the mixing of states and the corresponding 
energy levels and transition probabilities were obtained 
analytically in Refs. I10I14I As illustrated in Fig. [^a), 
there are a total of 2{2Ih + 1) quantum states (e.g. 20 
levels for bismuth with =9/2, 8 for arsenic which has 
Ih = 3/2). At high magnetic fields, the Zeeman states 
|z) = ms = ±1/2,= -Ih,-Ih±l, 

provide good quantum numbers. At lower fields, a new 
adiabatic set of states \i) = |±, m) must be employed, 
since ms and mj^ are not good quantum numbers, but 
m = ms ± mj,^ is. The relation between the Zeeman 
basis and the adiabatic basis is given by: 

|+,to) = cos ^ m — I) -I- sin ^ | —4,m -I- I) 

| — ,to) = — sin^|i,m — 4)-|- cos ^ | —m -I- |) , 

( 7 ) 

for all states except the two states with \m\ = Ih ± S 
(i.e. states 10 and 20 in Fig.j^a)) which remain unmixed 
at all fields. All other states mix with one other, form¬ 
ing doublets of constant m. The transformation between 
the Zeeman basis and adiabatic basis is given by simple 
rotation matrices 'Ry{j3m) and Ry(/3m)- Defining pa¬ 
rameters Xm = Ih{Ih ±l) — m? -1-1/4 and Zm ~ m-l- 
the angle of rotation is /3m = tsm~^[Xk/Zk\. 

For our system of interest, the Zeeman energy of the 
central system wg dominates over typical system-bath 
couplings Ja- This motivates a pure dephasing model 
(i.e. keeping only terms which don’t depolarize the states 
of the central system) whereby the bath dynamics is gov¬ 
erned by effective Hamiltonians depending on the state 
of the central system: ± ^bath- The 

key parameter of interest is: 

Pi{B) = = cos/3„, (8) 

which is the expectation value of the electron spin z- 
projection; it is no longer fixed at either ms = 1/2 or 
ms = —1/2 as for an unmixed qubit, but is instead a 
strongly field-dependent quantity Pi{B) S [—1/2 : 1/2]. 

OWPs correspond to operating the qubit at particular 
B values where Pu{B) ~ Pi{B). They provide a highly 
effective method of mitigating decoherence as illustrated 
in Fig. [^d). Recently, their importance has been recog¬ 
nised for certain donors in silicon such as arsenic or bis¬ 
muth. They were investigated theor etically in Refs. ITUF 
I12I14I and also experimentally,^^^^ extending the elec¬ 
tronic spin coherence time from 0.5 ms to 100 ms for an 
ensemble of donors in natural silicon. OWPs ha ve be en 
extensively investigated for classical field noise 
and in Ref. [55| dynamical decoupling was also studied. 


(a) OWP: = 4 G (b) :^OWP: 6 = 3200 G 




FIG. 2: (color online) Shows quantum many-body calcula¬ 
tions of the Hahn echo {N = 1) using the cluster correlation 
expansion (CCE) method, (a) Near OWPs, calculations us¬ 
ing a bath of independent spin pairs only (red, CCE2) do not 
even predict a finite decay time but, surprisingly, calculations 
with clusters of three spins (blue, CCE3) are already well- 
converged. The dashed lines used Eq. 0, a closed-form equa¬ 
tion derived from the short time behavior, found in Ref. E] 
to yield good agreement with experiments; this indicates that 
three-cluster results too give good agreement with measure¬ 
ments. Higher order CCE can encounter numerical diver¬ 
gences (which can be attenuated by ensemble averaging); this 
accounts for the discrepancies with CCE5. (b) Far from the 
OWP, independent pairs (CCE2) already give results in good 
agreement with CCE3-5 as well as experiments. The free in¬ 
duction decay (FID) is also shown for comparison. Note that 
Eq. (§ approximates the decay by a pure Gaussian. CCE 
calculations were performed for a bismuth donor in natural 
silicon for B along [100] and the |14) —>■ |7) transition for 
which Bowp = 799 G. In (a), B — 795 G while for (b), 
B = 3200 G. 

but Refs. Em considered quantum baths consisting of 
nuclear spins. 

The main eight OWPs for the bismuth system, occur¬ 
ring for R < 0.2 T are associated with avoided cross¬ 
ings. The transitions for four of these are shown by the 
arrows in Fig. Ha), while the other four correspond to 
forbidden transitions close by. The OWP we consider in 
this work is for the jl4) —)■ j7) transition which occurs 
at magnetic field Hqwp = 799 G. In practical realisa¬ 
tions, OWPs have become closely associated with field 
values where the derivative of the transition frequency 
/ = {Eu — El)/2Tr with respect to magnetic field vanishes. 
These are refereed to as “clock transitions” and corre¬ 
spond to suppression of classical field fluctuations, such 
as those arising from instrumental noise.^SlThe OWPs we 
consider are associated with reduced decoherence from 
quantum spin baths. Details of the difference between 
OWPs and clock transitions are given in Appendix [A] 

III. RESULTS: LOW-ORDER DYNAMICAL 
DECOUPLING 

A. Hahn spin echo (CPMGl) 

Understanding of decoherence for such mixed sys¬ 
tems nevertheless remains incomplete. In Refs. I1II12I591 
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CCE2 calculations were carried out to obtain coher¬ 
ence times for allowed and forbidden electron spin reso¬ 
nance (ESR) transitions coupling different pairs of states 
|u) —1^). These CCE2 calculations gave excellent agree¬ 
ment with experiment over most regimes. However, in 
the vicinity of the OWPs (where Pu — Pi), the CCE2 
Hahn echo decay failed to converge and no decay was 
obtained other than initially, for a short time. Single¬ 
central spin free induction decay (FID), in contrast, gave 
finite decays at all magnetic fields. 

In Ref. [12 an analytical expression estimating T 2 as 
a function of B was obtained, by inspection of the short 
time behavior of the form of the FID decays (which can 
be given analytically for each pair cluster): 

(9) 

The magnetic field dependence is wholly contained in 
the envelope, while the prefactor C{9) depends 

only on magnetic field orientation, the density of nuclear 
spins and their gyromagnetic ratio, but is independent of 
the strength of 

This simple closed-form equation gave remarkable and 
accurate quantitative agreement with experiment in all 
regimes, spanning orders of magnitude changes in T 2 , 
whether in the unmixed limit of a spin-1/2, or for cer¬ 
tain transitions which are ESR forbidden at high fields, 
or at OWPs. The universal validity of Eq. is wor¬ 
thy of discussion. Farther than about 100 G from the 
OWP, and where CCE is converged at CCE2, there is lit¬ 
tle difference between single-spin FID and Hahn echo de¬ 
cays; thus, it is not surprising that an equation obtained 
by considering the independent spin-pair contribution to 
FID can accurately model the Hahn echo experiments. 
Its validity within the OWP regions, however, is not yet 
fully understood. In particular, it remains unclear why a 
single C[9) prefactor suffices to accurately estimate ex¬ 
perimental T 2 , whether very far or very close to OWPs; 
and to describe different OWP regions (of which there 
are 16 for the bismuth system, with P^^i values varying 
by close to an order of magnitude). 

In the present work however, Eq. ([^ is employed 
simply to simulate the expected experimental behav¬ 
ior. In Fig. [2 we present converged CCE calculations 
for the Hahn spin echo (CPMGl) near the OWP. We 
show that including three-spin clusters (CCE3) give con¬ 
verged results while qubit-bath correlations from only 
spin pairs (CCE2) give little decay (red line) except at 
short timescales. However, it can also be seen that all 
orders have similar short time behavior and that the in¬ 
clusion of the three-clusters in effect recovers the short 
time behavior of the pair decays. Hence, the converged 
CCE agrees well with Eq. Q, which was derived from the 
early time decay of correlations from pairs and hence ac¬ 
counts for the experimental behavior observed in Ref [131 
This is one of the key results of this work. Details of the 
CCE calculations are given in Appendix [B} 


(a) OWP:e^„,-e = 4G 




Time (ms) 


FIG. 3: (color online) Shows dependence of the coherence on 
the number of dynamical decoupling pulses N, near (a) an 
optimal working point (OWP) and (b) far from an OWP, for 
modest numbers of N. (a) For B close to Bqwp, the T 2 times 
show comparatively little response to dynamical decoupling. 
Further, even though the initial coherence is extended with in¬ 
creasing N, the decays become ever more oscillatory. For low 
N, the independent pairs contribution is largely eliminated. 
Inset of (a): Showing complete suppression of the indepen¬ 
dent pairs contribution near an OWP; but showing also its 
gradual revival as N increases, (b) In contrast, far from the 
OWP, substantial (order of magnitude) enhancement of the 
T 2 time by dynamical decoupling is achieved with a moderate 
(preferably even) number of pulses. Decays for independent 
pair contributions (dashed lines, CCE2) and the converged 
quantum many-body numerics (solid lines, CCE4) are also 
compared, indicating that as A > 10, once again, the inde¬ 
pendent pair contribution is sufficient. Parameters are as in 
Fig.i and for CPMGA sequences. The converged GGE in 
(a) corresponds to CGE3. 


B. CPMGA with few pulses A ~ 2 — 20 

In Fig. we present comparisons of the response to 
dynamical decoupling near and away from an OWP by 
means of converged CCE calculations in both cases for 
N up to 16. One notable feature of the comparison is the 
insensitivity of OWP behavior to low numbers of pulses, 
in sharp contrast to the y^OWP regime where there is 
a factor of 3 “jump” in T 2 from CPMGl to CPMG2; 
while for larger A, we find T 2 ^ A as seen in Fig. [^c). 
However, OWPs are extremely effective at suppressing 
decoherence: for the point shown near the OWP, T 2 ~ 
100 ms already at CPMGl, while away from the OWP, to 
obtain comparable values, A ~ 100 pulses are required 
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as shown in Section V. 

Previous studies, including a recent study of the 
ESR dynamics of a phosphorus donor at X-band fre¬ 
quencies (a wstem without OWPs for electron qubit 
decoherence) observed a sharp increase in the coher¬ 
ence time between CPMGl and CPMG2. The spin 
pair contributions were also suppressed, requiring many- 
body correlations for convergence and thus exposing the 
latterP Nevertheless, GGE2 was shown to still give a rea¬ 
sonable approximation to the magnitude of the observed 
T 2 time, for both GPMGl and GPMG2. In the case away 
from an OWP, the FID is very similar to GPMGl. This 
is in contrast to the OWP, where GGE2 gives no decay 
at all, while the FID gave decay curves comparable to 
converged GGE3 (and Eq. ©)■ Thus, there is a drastic 
change from FID to GPMGl at OWPs; in contrast, for 
regimes away from an OWP, there is little change be¬ 
tween FID and GPMGl, but a strong enhancement for 
GPMG7V with > 1. 

The quantum numerics do evidence a clear dependence 
of the pair contribution on pulse number N. For exam¬ 
ple, in the inset of Fig. [^a), we have shown that, for 
a given field B in the vicinity of the OWP, as N in¬ 
creases to iV ~ 16, the pair contribution once again gives 
significant decay. To suppress decay for iV = 16 one 
must choose a value of B even closer to the OWP. In 
fact this is one of the main hndings of the present work: 
whether at OWPs or far from OWPs, our comparisons 
between many-body GGE3-5 and calculations involving 
only pairs show that increasing N gradually restores the 
importance of the pair contribution, relative to = 1 
OT N = 2, where many-body effects are seen to make 
the dominant contribution. In that case, and whenever 
independent pairs are dominant, we can employ well es¬ 
tablished, two-state analy tical pseudospin models of the 
qubit-bath dynamics. 

IV. ANALYSIS: PSEUDOSPIN MODELS 

We now proceed to analyse correlations from inde¬ 
pendent pairs in order to obtain insight on the effect 
of dynamical decoupling near and far from OWPs. As 
described above, assuming pure dephasing justified by 
Wo S> Ja, for the case of a pair of bath spins, the joint 
system-bath dynamics reduces to a simplihed two-state 
form for each of the two qubit states (upper and lower) 
and is governed by effective Hamiltonians: 

huj. — l^z “t” Ci2^x^-) (1^) 

where is a system-dependent detuning and 012 is the 
strength of the intrabath secular dipolar interaction. The 
Pauli matrices are written in the non-interacting bath ba¬ 
sis {lit) j ltt)}j where t (t) denotes spin up (down) for 
each of the two bath spins and we have ignored the trivial 
states Itt) and |tt) which cannot flip-flop. The detun¬ 
ings in our case are = Pu,i{Ji — J 2 ), where Ji and 


J 2 are the hyperfine coupling strengths for the two bath 
spins. The system-bath evolution is then obtainable ana¬ 
lytically from the eigenvalues uju,i = |-\/(A))()p^f 7 ^tL 2 p 
and eigenvectors of the above effective Hamiltonians. 

After preparing the initial qubit superposition, the 
GPMGV pulse sequence can be summarized as {T{t) — 
tt — T{t))^, with final evolution time t = Nt/2. The uni- 
taries T{t) represents free evolution and tt denotes the 
refocusing pulse which flips between |u) and |/): |m) (Z| -I- 
\l) (m| but leaves all other central states unperturbed. In 
the so-called pair-correlation approximation,the coher¬ 
ence decay is simply given by C{t) = J([^, where 

is the decay contribution from the fc-th spin pair 
and the product is over all spin pairs in the bath. 

For GPMGV, if 

(II) 

we can write for the Hahn spin echo case (i.e. GPMGl): 

= Aoi - iAu,i ■ d-, (12) 

where A„ = {A^, Ay, A^) and a is the vector of Pauli 
matrices in the bath basis: {lit) j Itt)}- The Au,i com¬ 
ponents depend on time and can easily be given explic¬ 
itly in terms of the pseudospin parameters; for exam¬ 
ple: AQ{t) = cos LUut cos u!it — sinu!ut sinujit cos {Ou — Oi) ; 
and Ay{t) = — sinsinw/f sin (0„ — 0/). where 0u,i = 
tan“^ (C 12 /Au,i). In fact, the only term which is not in¬ 
variant w.r.t. u I is Ay and thus A/ = (A^, —Ay, A^). 

The coherences oc |(H(0)|T)^ '^^Ti^^|H(0))| are 

obtained simply from L^^^{t) = 

For both GPMGl and GPMG2, the unitarity of the 
evolution of upper relative to lower states is broken by a 
term ex Ay. For GPMGl, 

LW(f) = l-2zAydyfW. (13) 

We can consider higher sequences; since 
and we obtain for GPMG2 

^ i _ 4iAy(A,d, - Axa,)fi^K (14) 

Both the above general expressions apply equally to ei¬ 
ther OWP or the ^OWP regimes. The only important 
difference between these regimes is that 9i for the 

approach to an OWP and 9^ = tt — 9i for the spin away 
from the OWP. Alternatively, from the explicit expres¬ 
sions for the components of A^.i, we see that the OWP 
condition is Ay —> 0; since Ay is the prefactor to both 
the above expressions, GPMGl and GPMG2 are equally 
suppressed at OWPs. 

For the thermal initial bath states |{,t) or |{,t)) the tem¬ 
poral coherence decays for the A:-th spin pair of the bath is 

i^r(t)i = I ati i;t) I = I (ni LfHt) in) i- The 

full decay is given by 

/:(f) = ni4'^^(t)l- 

k 


( 15 ) 
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We can easily obtain the coherence decay envelopes for 
CPMGl in general, assuming pulse interval r: 

= 2 t)\^ = 1 - AAlAl (16) 

where we drop the k label for convenience and = 
Ao(t), Ay = Aylr). For arbitrary even numbers of 
pulses, CPMGA^ such that N/2 is an integer. 


= 2Nt) = 1 



sin^ 


N(j){T) 

2 


(17) 


where cos(/)(r) = Aq(2t). An equivalent expression 
was obtained in Ref. |6l Both expressions Eq. ( 16 ) and 
Eq. ( [l7| are equally valid for both regimes (OWP and 


A. OWP limit 

The only important difference between these regimes is 
that —>■ 9i for the approach to an OWP and 6 ^ = tt —9 1 
for the spin away from the OWP. Alternatively, from the 
explicit expressions for the components of we see 
that the OWP condition is \Ay\ 0. Thus, the sup¬ 
pression of qubit-bath correlations from pairs for OWPs 
is of the same order for GPMGl, GPMG2 or any other 
even-pulsed GPMG: for all bath spin pairs equally, the 
decay due to correlations from an independent pair uni¬ 
formly tends to zero as {Ay^kY 0 f^e fc-th pair as 
B —>■ Bowp- 

The dependence on N is entirely contained in the 
s\T? N / 2 term. If N(f>{T) 1 then increasing N 

has a strong amplifying effect on the signal, while if 
N(I){t) ^ 1, increasing N simply results in oscillatory 
behavior. Near OWPs, from the expression for Ao(2r), 
we see that if 9^ = 9i, 4>{T)f2 ~ (a;„ -I- uji)t. Hence we 
only expect a response to dynamical decoupling if r is 
sufficiently small (i.e. if r < (a;„ -t-w;)”^). 


B. /OWP limit 

In contrast, for GPMG away from an OWP, the A^ 
prefactor is still there, but is not small. The origin of 
the suppression of correlations from independent pairs 
for small numbers of pulses is more subtle to analyse 
with the pseudospin model. For GPMG2 (/OWP limit), 
we obtain: 


|/:(2)(t)|2 = i_64A2/2^t (18) 

The large jump in T 2 from GPMGl to GPMG2 was also 
analysed in Ref. [HI In the notation of Ref. [HI we see that 
for GPMGl, the decay envelope is of order n^, while for 
GPMG2 it is of order n^nl, where Ux = sin0^l = sin 61 
while riz = cos 9u = — cos 9i. Since the bath spans all 
angles \9u,i\ = [0,7r/2] one cannot a priori assume sin0^i; 
is small. However, previous numerical studies support 


(a) OWP:e„„p-e=4G 




Time (ms) 

FIG. 4: (color online) Shows behaviors of coherence decays 
for large numbers {N) of dynamical decoupling pulses (a) 
near and (b) far from OWPs; as shown in Fig. |^a), in this 
regime, correlations from independent pairs once again dom¬ 
inate the decays in all regimes so CCE2 is converged and 
plotted. The behavior at OWPs is now sensitive to N but the 
decays here become increasingly oscillatory as N and T 2 both 
become large; we attribute this to large numbers of bath spin- 
pair frequencies becoming resonant with the pulse spacing. It 
indicates the behavior one might expect in a single-shot sin¬ 
gle spin study. The smooth lines are fits to the decays and 
indicate the behavior expected after ensemble averaging. Pa¬ 
rameters are as in Fig. 


the idea that those spin pairs which have |Ji — J 2 I ^ 
\Ci 2 \ (i.e. are strongly coupled to the central system) and 
therefore small pseudospin angles, dominate the Hahn 
echo contribution.li^ For GPMG2, such strong-coupled 
spin pairs are strongly suppressed, and so T 2 becomes 
dominated by more weakly coupled spin pairs which are 
less effective in decohering the qubit. 


V. RESULTS: CPMGA WITH LARGE N 

For large N, decays from independent pairs only 
(GGE2) are restored as well as the sensitivity to dynam¬ 
ical decoupling at OWPs. Even for A = 16 (Fig. 
we see that the initial period of no decay C{t) ~ I is 
prolonged. For larger N (Fig. |^, the enhancement of co¬ 
herence even at OWPs is clear, but however, the decays 
become extremely noisy. The noise can be attributed to 
the timescales of individual nearby spin clusters and the 
time interval between pulses. For these long coherence 
times (^ Is) there are very large numbers of resonances. 




































The CPMG sequence provides a means of amplifying 
noise from nearby clusters whenever pulse intervals be¬ 
come resonant with the characteristic cluster frequency. 
While this makes the CPMG a valuable technique for 
spin detection!^ large numbers of such resonances are 
undesirable if the aim is to protect qubit coherence. In 
contrast, far from OWPs, the decays for high N remain 
relatively smooth. While the noise at OWPs can be mit¬ 
igated by ensemble averaging, this is likely introduce a 
considerable disadvantage in terms of single-shot opera¬ 
tion of a single qubit. 

VI. DISCUSSION 

In sum, we have seen that a key difference between 
OWP and ^OWP behaviors arises from the Ay oc 
sin^(6i„ — 6 i) prefactor which globally suppresses all inde¬ 
pendent pair contributions on the approach to an OWP, 
and accounts for the drastic effect at OWPs, but which 
is independent of N and has little effect far from OWPs. 
However, to analyse decays resulting from dynamical de¬ 
coupling one must consider the remainder of the expres¬ 
sion in Eq. 0. which reflects the dependence on N. 
The ineffectiveness of dynamical decoupling near OWPs 
for small N can also be understood with an intuitive pic¬ 
ture considering the relevant timescales of the system. 
For dynamical decoupling to be effective, the time inter¬ 
val between pulses {t/2N) must be shorter than to the 
correlation time of the bath Tc- Since typical intra-bath 
interactions are at most a few kHz, Tc ~ 1 ms. Near 
the OWP, ~ w* and ~ 0;, so the frequency of the 
bath noise spectrum uiu,i) is appreciably higher than 
1 /tc and thus dynamical decoupling becomes ineffective 
in extending the coherence time T 2 'A> Tc- At short times 
and for high N however {t/2N < tc), dynamical decou¬ 
pling does protect the central system as evidenced for 
CPMG16 in Fig. |^a) and higher N in Fig. In con¬ 
trast, dynamical decoupling is far more effective in ex¬ 
tending T 2 away from the OWP and for relatively small 
N (Fig. [^b)); although the pseudospin frequencies are 
comparable, the pseudospin fields are In opposing direc¬ 
tions (0„ ~ TT — 0;), thus, the frequency of noise is much 
slower and becomes comparable to 1 /tc ~ 1 /T 2 . 

Finally, it is important to note that for direct quan¬ 
titative comparisons between our dynamical decoupling 
calculations and experimental ensemble measurements. 
Inhomogeneous broadening due to ^®S1 nuclei might also 
have to be factored in (see Appendix . 

VII. CONCLUSIONS 

Understanding the interplay between optimal working 
points and dynamical decoupling Involves understanding 
of the quantum behavior as a function of the two limits 
B —>■ Bqwp and V —)■ 00 corresponding to approach¬ 
ing an OWP and simultaneously Increasing the number 


of dynamical decoupling pulses. An underlying question 
of physical interest is when decoherence is the result of 
the magnetic noise from independently flip-flopping pairs 
of spins and when consideration of the many-body na¬ 
ture of the quantum bath is important. The answer is 
of practical importance. For one, if decoherence is due 
to flip-flopping pairs, there are widely used models (such 
as the analytical pseudospin expressions in Section IV) 
which can be used to accurately calculate decays. Oth¬ 
erwise, more complex full many-body numerics become 
essential to simulate and fully understand experimental 
behaviors. The clear answer Is that for low order dynam¬ 
ical decoupling, the elimination of correlations from in¬ 
dependent pairs is so drastic at OWPs, that many-body 
numerics is almost Indispensable for full understanding 
and accuracy. Even away from OWPs, it was shown that 
many-body effects make a large contribution for N < 10. 
However, once V —> 00 we find little difference between 
Independent-pair and many-body results. 

For practical applications, by solving for the many- 
body qublt-bath dynamics to calculate coherence times, 
one can hope to Identify the best strategy for enhanc¬ 
ing the coherence of donor qubits whilst still keeping the 
nuclear spin bath of naturally occurring silicon for its 
potential technological use. By operating near OWPs 
without dynamical decoupling, the maximum achievable 
T 2 is 0.1 s due to in homo geneous broadening from the 
environmental nuclei .1^21131 Pqj. isotopically purified sam¬ 
ples in which the nuclear spin bath is nearly eliminated, 
T 2 at the OWP was measured to be about 1 s and is lim¬ 
ited by decoherence mechanisms involving donor-donor 
interactions.^ Therefore, to bridge this single order of 
magnitude difference in T 2 at OWPs without resorting 
to isotopic purification, dynamical decoupling should be 
applied with at least a few hundred pulses. The effect of 
dynamical decoupling in extending coherence times near 
an OWP is marginal with a moderate number of pulses 
(up to V ^ 16) In contrast to the usual regimes far from 
OWPs. For high donor concentrations, the timescale 
of donor-donor decoherence is comparable to the T 2 ob¬ 
tained In a nuclear spin bath, hence one might also want 
Investigate suppressing those mechanisms with dynami¬ 
cal decoupling. 

However, combining dynamical decoupling with OWPs 
is not without its drawbacks. As T 2 ,N —>■ 00 , poten¬ 
tially many spins in a silicon bath may become reso¬ 
nant with the dynamical decoupling pulse spacing, re¬ 
sulting in very noisy decays in single central spin realisa¬ 
tions. Although ensemble measurements are unaffected 
by this noise, this means that for single-qubit operations, 
if OWPs can be exploited, their extraordinary potential 
for coherence suppression may be sufficient. 
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Appendix A: OWPs vs. clock transitions 


The expansion is given by 


K 


^K{t) 


T\tI(ZK 


(Bl) 


The OWPs we consider correspond to suppression of 
decoherence in quantum spin environments. Similar 
points of in terest are clock transitions or df/dB = 0 
points,li^E^ where decoherence arising from classical field 
noise is suppressed. For donor spin systems, the OWP 
is close to but not exactly at the df/dB = 0 point, and 
not all df /dB = 0 points are OWPs as shown in Ref. [121 
This is because the donor is coupled to the nuclear spin 
bath primarily via the electron, and the couplings be¬ 
tween the host donor nucleus and bath nuclei are negligi¬ 
ble compared to the electron-bath hyperfine interaction. 
It is easy to show that the fields at which df/dB = 0 
satisfy 


0 — Pu — Pi + 


SAm 

TTj’ 


(Al) 


where 6 ~ 10“^ is the ratio of the electronic to host 
nuclear gyromagnetic ratios and Am is the difference in 
quantum number m between the upper and lower levels. 
Since <5 1, the OWP (where Pu = Pi) and df/dB = 0 

points nearly coincide. 

It is important to note that all of the transitions with 
OWPs couple two neighbouring avoided crossings. Selec¬ 
tion rules were detailed in Ref. dH but all such transi¬ 
tions have Am = ±1 which implies that {u\Sz\l) = 0. 
Thus, magnetic field fluctuations represent pure dephas¬ 
ing noise. One might also consider the possibility of cre¬ 
ating a superposition of two states lu) and |Z) at a single 
avoided crossing; for example, the superposition |II)-|-|9) 
in Fig. at the avoided crossing between these states 
at B = 0.2IT. Although the |II) —>■ |9) transition is never 
allowed, such a superposition might be created by a two 
pulse excitation from level |I0). Both states are at zero 
energy gradient {dEuy/dB = 0) so coherences are to first 
order insensitive to dephasing noise; however, as shown in 
Ref. [in in that case {u\Sz\l) 0 so magnetic fluctuations 
couple the states in the superposition and thus coherence 
is vulnerable to depolarisation by magnetic noise. 


Appendix B: Cluster correlation expansion 
simulations 

The CCE is a well-established method for accurately 
calculating the coherence decay £{t) of a central spin 
system in a quantum spin bath.^^ In the CCE formula¬ 
tion, the spin bath is decomposed into groups of spins or 
“clusters” and the closed evolutions of clusters interact¬ 
ing with the central system are combined to approximate 
the exact coherence involving the entire spin bath. 


where A is a subset of the bath and the first product is 
over all subsets of the bath. The irreducible or “true” cor¬ 
relation term for cluster 1C, C]c(t), is recursively defined 
as follows. First, the coherence is calculated by exactly 
solving for the combined qubit-bath dynamics governed 
by the total Hamiltonian in Eq. Q, but including only 
those bath spins contained in 1C. Second, the resulting 
Cjcit) is divided by all correlation terms formed out of 
the proper subsets TZ of JC. The usefulness of the CCE 
becomes clear when truncating the expansion, 

C[k]{t)= ll £dt), (B2) 

|K|<fe 

in which only clusters containing a maximum of k bath 
spins are included. The lowest non-trivial order for nu¬ 
clear spin diffusion is the pairwise correlation approxima¬ 
tion (fc = 2 or CCE2), given by Eq. ( |I^ , and we calculate 
up to The coherence decay is often considered 

converged when <C l,Vf and for our 

case this condition is satisfied for up to fc' = 3. 

Note that in the CCE method, the pure dephasing ap¬ 
proximation is not required, and the interaction Hamilto¬ 
nian in general includes terms which depolarize the states 
of the central system. Our CCE calculations include the 
S~ +S^ I~ terms in the hyperfine interaction Hamilto¬ 
nian Eq. (§ , but we find that these give small corrections 
to the case when only terms are included. This was 
expected due to the large mismatch between electronic 
and nuclear gyromagnetic ratios. 

For our simulations, crystal sites of a cubic silicon 
superlattice were uniformly populated with ^®Si nuclei 
(J = 1/2) with the natural fractional abundance of 0.0467 
and with equal probability of spin-up and spin-down for 
the initial states. To calculate the electron-bath hy¬ 
perfine couplings, we use the Kohn-Luttinger electronic 
wavefunction for the bismuth donor in silicon with an 
ionization energy of 0.069 eV.^^ The total size of the spin 
bath is dictated by the spatial extent of the wavefunc¬ 
tion which decays exponentially with distance from the 
donor site, and a superlattice of side length 160 A(with 
10^ impurities) gave convergent coherence decays. Due 
to cubic decay of the dipolar interaction as the distance 
between a pair of ^®Si spins is increased, it is not neces¬ 
sary to include all spin clusters in the calculation. At the 
lowest non-trivial CCE order (CCE2), spins separated 
by at most the 4-th nearest neighbor distance in silicon 
(•\/TTao/4, where oq = 5.43 A) gave convergent decays 
away from OWPs. To choose 3-clusters (i.e. including 
three bath spins), we loop over all sites in the crystal 
and add to each 2-cluster only those spins that are at 
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most separated by vTTao/4 from any of the two spins in 
the 2-cluster. The same procedure was applied to choose 
higher-order clusters, by adding spins to clusters one or¬ 
der down. In Fig. and Fig. the numbers of 2, 3, 4, 
and 5-clusters were each of order « 10"^. 



Magnetic field, B (G) 


FIG. 5: (color online) Sharp B-field dependence of T 2 for var¬ 
ious CPMG orders near an OWP. Inhomogeneous broadening 
from ^®Si nuclei can be incorporated by convolving the decays 
with a Gaussian B-field distribution centred about B (here 
centred about 797 G) and with standard deviation w ~ 2 G 
(dashed line). For a donor concentration of 3 x 10^® cm'®, T 2 
is limited by donor-donor processes at about 300 ms.^ The 
T 2 lines were calculated for bismuth donors using the GCE up 
to 3rd order and for B || [110]. The OWP under investigation 
is shown in red at 799 G. 


Appendix C: Effect of inhomogeneous broadening 

For direct quantitative comparisons between our dy¬ 
namical decoupling calculations and experimental ensem¬ 
ble measurements, inhomogeneous broadening might also 
have to be considered. As shown in Fig. for various or¬ 
ders of CPMG, T 2 varies sharply with magnetic field over 
a few G near an OWP. Inhomogeneous broadening of B 
due to ^®Si impurities has a FWHM of about 4 G in nat¬ 
ural silicon and may therefore need to be included in the 
calculation in order to predict the shape and rate of ex¬ 
perimental decays near OWPs. The broadening can be 
simulated by convolving the decays Csit) with a Gaus¬ 
sian magnetic field distribution with standard deviation 
w ~ 2 G: 

Dsit) = —= / e 2^2 CB{t)dB'. (Cl) 

wv27r J 

Depending on the donor concentration, donor-donor pro¬ 
cesses may also need to be included. For example, for a 
donor concentration of 3 x 10®® cm'®, T 2 near an OWP 
is limited by direct flip-flops of the central donor with 
other donors in the ensemble.l^ The measured T 2 in an 
isotopically purified sample ranges from 0.2 — 2 s (T 2 
for a donor concentration of 3 x 10®® cm'® is 300 ms). 
Therefore, care should be taken to include donor-donor 
processes very near the OWP (within about 1 G), where 
nuclear spin diffusion coherence times are comparable to 
those of donor-donor processes. 
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